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Abstract 

In this paper, we discuss the constraint on the relic gravitational waves by both temperature 
and polarization anisotropies power spectra of cosmic microwave background radiation. Taking 
into account the instrumental noises of Planck satellite, we calculate the signal-to-noise ratio S/N 
by the simulation and the analytic approximation methods. We find that, comparing with the 

o ' 

BB channel, the value of S/N is much improved in the case where all the power spectra, TT, 
TE, EE and BB, are considered. If the noise power spectra of Planck satellite increase for some 

\ reasons, the value of S/N in BB channel is much reduced. However, in the latter case where all 

the power spectra of cosmic microwave background radiation are considered, the value of S/N is 
less influenced. We also find that the free parameters A s , n s and nt have little influence on the 

. value of S/N in both cases. 

o 
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I. INTRODUCTION 



A stochastic background of the relic gravitational waves (RGWs), generated during the very 



a. 



curly slag* 1 of the I'nivcrse lj], is a necessity did aUx'l by general rclat ivil.y and quantum mechanics 



4, la] , and their detection plays a double 



The RGWs have a wide range spreading spectra 
role in relativity and cosmology. 

One of the important methods for the detection of RGWs is by the cosmic microwave background 
(CMB) power spectra, including the temperature anisotropies (TT) power spectrum, the polariza- 

Onnnn 

tion (EE and BB) power spectra, and the cross-correlation (TE) power spectrum [a. 1 71. la. 19. llOQ. 
By observing the lower-order CMB multipoles, one can detect the signal of RGWs at the very 
low frequencies (10 -17 ~ 10 -15 Hz). One has to note, besides the gravitational waves of quantum- 
mechanical origin 2j , the classical gravitational waves were generated at later stages of cosmological 
evolution However, their wavelengths are much shorter than the present Hubble radius and 
therefore they do not affect the lower-order CMB multipoles. 

As well known, the CMB has certain degree of polarization generated via Thompson scattering 
during the decoupling in the early Universe [121 ] . In particular, if the RGWs are present at the 
photon decoupling in the Universe, the magnetic type of polarization (i?-polarization) will be 
produced (j, 8, [9]. This would be a characteristic feature of RGWs on very large scale, since the 
density perturbations will not generate this polarization. So a natural way for the detection of 
RGWs is by observing the signal of U-polarization of CMB. This is the so-called "55" method. 
However, the amplitude of the U-polarization is expected to be very small. In addition, the 
jj-polarization is prone to degradation by various systematic effects on a wide range of scales 
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161 ]. The current 5-year Wilkinson Microwave Anisotropy Probe (WMAP5) observation 



only gives an upper limit £(£ + l)C£f 2 _ 6 /2TT < 0.15^K 2 (95%C.L.) 



such as the Planck 



lfiJl, Clover [19], Spider [20j], QUITE 



171 ] . The forthcoming projects, 



211 ] . are expected to be much more 



sensitive for the detection of the CMB 5-polarization. 

Due to the disadvantage of ll BB" method, it is necessary to look for the new method for the 
detection of RGWs in the CMB. In the previous work [9J, the authors found that, in the large scale 
(£ < 50), the RGWs generate the negative TE spectrum. However, if the TE spectrum is generated 
by density perturbations, it should be positive. This suggests that the signal of RGWs can also 
be detected by the CMB TE spectrum. Comparing with the 5-polarization, the amplitude of the 
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22, 



23J, the authors have developed several 



TE spectrum is nearly two order larger. In the works 
ways to detect the signal of RGWs directly from the CMB TE spectrum. These are the so-called 
"T£" method. 

In the work 23j we found the WMAP5 TE data contains a hint of the presence of RGWs 
contribution. In terms of quadrupole ratio R, the best-fit model produced R = 0.24, which cor- 
responds to the tensor-to-scalar ratio r ~ 0.48. Because of large residual noise, the uncertainty 
of this determination is still large. We also found, if considering the Planck instrumental noises, 
U TE" method can detect the signal of RGWs at 2a level when r > 0.3. If considering the ideal 
case with full sky and no noise, a TE 7 method can detect the RGWs at 2a level when r > 0.1. By 
comparing the detection abilities of 11 TE 7 and ll BB" methods, we found that, taking into account 
of the instrumental noises of Planck satellite, ll BB" method is more sensitive when r is small. 
However, if the noise power spectra or the amplitudes of RGWs increase, the sensitivity of " TE' 
method becomes better than that of a BB" method. 

I — I 

In this paper, we shall extend the 11 TE 7 method in our previous work [23J]. By calculating the 
values of S/N, we shall compare the detection abilities for the RGWs in the following four cases. 
The first case ("fl" case) is the so-called il BB" method, where only the S-polarization spectrum 
is considered. In the second case, we include not only the CMB TE spectrum, but also the TT 
spectrum. We call it as the "CT" case ("C" standing for the cross-correlation power spectrum 
and "T" standing for the temperature anisotropies power spectrum). We shall expect that, taking 
into account the contribution of TT spectrum, the detection ability will be much improved. In the 
third case, besides TE and TT spectra, we also include the EE spectrum. We call it as the "CT-E" 
case. The fourth one is the "CTEB 7 case, where the contributions of TE, TT, EE and BB power 
spectra are all considered. The detection ability in this case is expected to be much more sensitive 
than the other three cases. 

The organization of this paper is as follows. In Section II, the primordial power spectra of 
RGWs and density perturbations, the CMB power spectra and the corresponding estimators are 
introduced. In this section, the probability density functions (pdfs) for the estimators are also 
discussed. In Section III, we introduce four ("5", 'CT", "CT£" and U CTEB') cases for the 
detection of RGWs in the CMB. The likelihood functions are also given this section. In Section IV, 
by constructing the likelihood functions based on the simulated data, we shall investigate the values 
of the signal-to-noise ratio S/N in these four cases. This quantity describes the detection abilities 



of RGWs in the different cases. We firstly introduce the simulation method. In the simulation, 
when constructing the likelihood functions, we only consider one free cosmic parameter, the tensor- 
to-scalar ratio r. We find when r > 0.06 (r > 0.16, r > 0.13, r > 0.05), the signal of RGWs can be 
detected in "B" ("CT', "CT£", "CTEB") case at 2a level. In Section V, we discuss the analytic 
approximation of the likelihood functions. By the analytic approximation formulae, we obtain a 
simple analytic form of S/N, which clearly shows the dependence of S/N on the amplitude of 
RGWs and the noises. By analyzing the analytic form of S/N, we find that, in "CT" and a CTE" 
cases, the main contributions come from the data in the intermedial scale 10 < i < 70. However, 
in U B" case, the main contribution comes from the data in the very large scale t < 10. In Section 
VI, by the simulation method, we find that, the free parameters (the tensor spectral index i%t, 
the amplitude of scalar spectrum A s and the scalar spectral index n s ) have little influence on the 
determination of RGWs. Section VII is the conclusion that summarizes the main results of this 
paper. 



II. GRAVITATIONAL FIELD PERTURBATIONS, CMB POWER SPECTRA 
AND THEIR ESTIMATORS 

A. Primordial power spectra of the gravitational field perturbations 

The CMB temperature and polarization anisotropies power spectra are determined by the pri- 
mordial power spectra of density perturbations (scalar perturbations) and RGWs (tensor perturba- 
tions), and the time evolution of these perturbations during and after the epoch of recombination. 
Before proceeding with the CMB power spectra, it is necessary to introduce the primordial pertur- 



bation spectra, which are usually assumed to 



)e power-law. This form is a generic prediction of a 



wide range of scenarios of the early universe 



24]. In general there might be deviations from a 



25J]), but 



power-law, parametrized in terms of the running of the spectral index (see for example 
we shall not consider this probability in the current paper. Thus the primordial power spectra of 
the perturbation fields have the forms 

p.(k) = Mk ) (£) , p t (k) = am (£) , (i) 

where n s and nt are the scalar and tensor spectral indices, respectively, kg is pivot wavenumber, 
which can be arbitrarily chosen. In the WMAP paper [261 ] . the pivot wavenumber fco = 0.002Mpc _1 



is used, which is close to the observable horizon. The scale ko = 0.05Mpc _1 is also commonly used, 
being the default scale of CAMB package j^]. A number of authors have suggested other pivot 
wavenumber k$ for different reasons [28]. In Eq. ([T]), A s (ko) and At(ko) are the amplitudes of the 
primordial scalar and tensor spectra respectively, at the pivot scale k$. 

We can re-parameterize the tensor power spectrum amplitude At(ko) by the "tensor-to-scalar 
ratio" r, which is defined by 

ar quadrupole R is also quoted when refer- 



23 



29j). The relation between R and r is 



In addition, the ratio of tensor quadrupole to seal; 
ring to the tensor-to-scalar ratio (see for instant 
somewhat cosmology-dependence, especially on the dark energy density The conversion is 
R ~ ° § i 04-0 82^A+2£^ nA r fity - For the cosmological models with = 0.76, these two definitions 
are simply related by r ~ 2R. In the previous work [23|], we have adopted R. However, in this 



paper we shall use r, the default quantity used in the CAMB package 



Using Eqs. (P) and ([2]), one can evaluate r at a different wavenumber k\, 

r(h) = r(k ) . (3) 

In the following discussion, we shall discuss the constraint on the parameter r by the simulated 
data. From the relation ([3|) we find that, if the spectral indices nt and n s are fixed as = 
and n s = 1 in the likelihood analysis (the case in Sections IIVI and |V|) . we have r{k\) = r{k$), the 
tensor-to-scalar ratio r has the exactly same value at all pivot wavenumbers. So we do not need to 
differentiate the tensor-to-scalar ratio at the different pivot scales. However, if the spectral indices 
are free parameters in the likelihood analysis (the case in Section |VT|) . comparing with r(fco), the 
constraint of r{k\) is also influenced by the spectral indices. So in this case, we should differentiate 
the tensor-to-scalar ratio at the different pivot wavenumbers. This effect will be clearly shown in 
Section EH 

We should mention that, in the specific early universe models, the parameters A s , n s , r and nt 



are always not separate 



321 ] . However, in this paper we shall avoid any specific model, and consider 



the parameters A s , n s , r and nt as the independent parameters in the data analysis. 
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B. CMB power spectra 



The CMB radiation field is usually characterized by four Stokes parameters (I ,Q,U,V). I is 
the total intensity of radiation, Q and U describe the magnitude and direction of linear polar- 
ization, and V is the circular polarization. From these Stokes parameters, we can construct four 
invariant quantities (I ,V, E, B), which can be expanded over ordinary spherical harmonics (see 
[j| for details). The set of multipole coefficients (af m , aj m , af m , af m ) completely characterize the 
intensity and polarization of the radiation field. Since Thompson scattering of initial unpolarized 
light cannot generate circular polarization, we shall not consider the V stokes parameter in the 
following discussion. 

In general, the output of the CMB experiment a x m (X = T, E, B), consists of two contributions, 
the signal convolved with the beam window function af m (s)We and the noise a«L(n), i.e. 

a? m = a i m ( s We + af m (n). (4) 

We shall use the notations (s) and (n) to denote the signal and noise. These two contributions are 
uncorrelated to each other. 

Assuming the primordial perturbation fields (including the scalar perturbations and tensor 
perturbations) are Gaussian fields, which induces that the signal term af m {s) has the covariance 



23] 



\{af m { s ) a fm'( s ) + a fm( s )4l'(*)) = C xx '8 u ,8 mm ,, (5) 

where C xx ' is known as the CMB power spectra, which depends on the cosmological inputs. 
When X = X', C xx ' is the auto-correlation power spectra, and when X ^ X' , C xx ' is the 
cross-correlation power spectra. In absence of any parity-violating processes, the only survived 
cross-correlation power spectra is CT E 



33J. So the temperature and polarization anisotropies can 
be described completely by four power spectra: Cf T , Cf E , Cf B and Cj E '. 

The noise terms CLf m {n) and the window function We depend on the experiment. We assume 
the noise is a spatially uniform Gaussian white noise. For an experiment with some known beam 
width and sensitivity, the noise power spectra and window function can be approximated as 

N l X ' = \{ a f m ( n ) a &m'{ n ) + a ¥m( n ) a *L'(. n )) = (#FWHM<7x) 2 <5xX'<5«'<5mm', (6a) 



Wi = exp 



1(1 + 1) flp WHM 
2 81n2 



(6b) 



where #fwhm is the full width at half maximum of the Gaussian beam, and ax is the root mean 
square of the instrumental noise. Non-diagonal noise terms (i.e., X ^ X') are expected to vanish 
since the noises contributions from different maps are uncorrelated. The assumption of a spatially 
uniform Gaussian noises spectrum ensures that the noise term is diagonal in the I basis. In this 
paper, we shall consider the Planck instrumental noises. There are several frequency channels for 
the detection of CMB in Planck satellite 



181 ] . In this paper, in order to estimate the Planck noises, 
we only adopt the frequency channel at 143GHz, which has the low foreground levels and the lowest 
noises power spectra. In this channel, we have 

#fwhm = 7.1', &t = 6.0/xK, oe = &b = 11.5//K. (6c) 

Inserting these into Eq. ()6a|) . we obtain the noise power spectra 

Nj T = 1.53 x KTVk: 2 , Nf E = Nf B = 5.58 x lO" 4 ^ 2 , Nj E = 0. (6d) 

Considering Eqs. fl3J), and (f6a|) . we obtain the covariances of the terms af m , which are 

^^trnP-Vm' + a fm a Vm') = ( C f X Wf + N XX )5 U /5 rnm i. (7) 

C. Estimators of the CMB power spectra 

In Section III Bt we have introduced the CMB power spectra, which are defined as ensemble 
averages over all possible realization of the CMB field. However, in CMB observations, we only 
have access to one single realization of this ensemble. In order to obtain information on the power 
spectra from a single realization, it is desirable to introduce the estimators of the power spectra, 
which are observable quantities. 

In the full sky case, and taking into account the noises, the best unbiased estimators Df x ' for 
the CMB power spectra Cf- X ' are defined by 



3, 



^xx 1 _ ( 1 ST (n x x>* , x*x> 
\ m=-£ 



a x A " i 1 T, («;<-/'?»'/ + • ' ■ £) - 'J wf , (8) 

where n is the number of the degree of freedom for a fix multipole I. In the full sky case, we have 
n = (21 + 1). The expectation values and the standard deviations of these estimators are 



XX' _J(Cf x +N^W e - 2 )(Cf x '+Nf' x 'W e - 2 )HCi cx '+N^'W e - 2 ) 



(9) 
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It is necessary to investigate the pelfs for D xx ' , which have been derived in 23], based on the 
assumption: the primordial perturbation fields and noise fields are independent Gaussian fields. In 
this subsection, we shall briefly introduce the results as follows (the similar results are also obtained 



in the Refs. [35l . 



m 



The pdf of the auto-correlation estimator D¥ x is known as the x 2 distribution, which is 



f(D xx ) = W 2 )^- 2)/2 e- V/2 (10) 
> 2^T{n/2){C xx Wf + N xx Y 1 ' 

where n = (21 + 1) is the degree of freedom for the multipole 0. in the full sky case. The quantity 

V is defined by V = n{D xx Wf + N XX )/{C XX W 2 + N xx ). 

The joint pdf for the estimators Dj E , Dj T and D EE is the following Wishart distribution 

f (n TE n TT n EE } _( 1 ) W/2 (nW?nxy-z*)W> 

n * ' 1 ' 1 > ~ \ A{1-Pl)(ajafy] 7rV2r(n/2)r((n-l)/2) (n) 

x ex P { " 2(T i ^y ( J^f + J^fy ~ 3^e) } , 
where the quantities x,y, z are defined by: x = n(Df T W 2 + Nf T ), y = n{Df E Wf + Nf E ), 
z = nDj E Wf. af = ^Cj T Wf + Nf T , of = \jc EE W} + N EE are the standard deviations of 
the multipole coefficients aj m and af m , respectively, pi is the correlation coefficient of aj m and af m , 
which can be written as, 

r<TE 

(12) 



\CJ T + Nf T W £ - 2 )(C EE + Nf E Wf) 

From the Wishart distribution (jlip . we can derive the joint pdf of the estimators Dj E and Dj T 
by integrating the variable Df E , the finial result is 

f(Dj E ,Dj T ) = (nWffx^ {21+^(1X1 - pi)(ajr(* E f}^ 

We can also obtain the joint pdf for all the four estimators: Dj E , Dj T , D EE , Df B . Since B- 
polarization estimator Df B is independent of the estimators Dj E , Dj T and Df E , the total joint 
pdf is the product of the Wishart distribution f(Dj E , Dj T , Df E ) in (jlip and the x 2 distribution 
f(Df B ) in ([10]) with XX = BB, i.e. 

f(Dj E , Dj T , D EE , Df B ) = f(Dj E , Dj T , D EE )f{Df B ). (14) 

We should notice that, the above results are all based on the assumption of full sky coverage. 
However, real experiments can only see a fraction of sky. Even for satellite experiments, a map cut 



S 



must be performed in order to eliminate point sources and galactic plane foreground contamina- 
tions. As a result, different multipole moments af m become correlated with each other The 
exact pdfs of the estimators in this case takes a rather complicated form, depending on the shape 
of remaining observed portion of sky [35f] . However, for experiments probing almost the full sky 
(e.g. COBE, WMAP, or Planck), correlations are expected only between neighboring multipoles. 
In order to simplify the problem, one can take a«L's to be uncorrelated, and introduce a factor 



23, 



for the estimators with 



/sky; which denotes the observed fraction of sky. As was shown in 
the multipole number £, the number of degree of freedom reduces to n c g = (21 + l)/ s ky (instead of 
n = 2£ + 1). Thus, compared to the full sky, the inclusion of cut sky reduces the degree of freedom 
in the definition of the estimators D¥ x ' . In this work, we shall discuss the CMB field with the cut 
sky factor 

/sky = 0.65, (15) 

which is suggested by Planck bluebook 



181 ] . In all the following discussion, we should remember 
to replace n with the effective degree of freedom n e g, when using the result in and the pdfs in 

(Hod, dm a or m 



III. FOUR CASES TO DETECT RGWS IN THE CMB 

In this paper, we shall investigate the detection abilities for the RGWs in the following four 
cases: "5" case, "CT" case, "CTE J case and "CTEB" case, which will be introduced separately 
in this section. 

A. "5" case 

The first case is the well-known "55" method. In this case, one can detect the signal of RGWs 
only by the observable D BB , which satisfies the x 2 distribution in Eq. (|10p . 

In order to study the determination of cosmic parameters from the observed data, we shall 
consider the likelihood function. The likelihood is a term, customarily, used to call the probability 
density function considered a function of an unknown parameter. Up to a constant, independent 
of its arguments, the likelihood is defined as the pdf of the set of the moments Df B given Cf B , 
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i.e. 



£ B *Hf(D? B ). (16) 

i 

Using the pdf in (fTUj) . and considering the effective degree of freedom n e ff = (2£ + l)/ s k y in the cut 



sky, the likelihood function in (|16p can be rewritten as 

f ( Df B + N, BB W,- 2 \ ( D? B + N? B Wf 2 \\ , , 

where Ci is the constant for the normalization. The noise power spectrum N BB for Planck mission 



is given by Eq. (|6d|) . 

Since the 5-polarization can only be generated by the gravitational waves, the observable 5- 
polarization power spectrum includes a clean information of the gravitational waves. This is the 
advantage of "55" method. However, the amplitude of the 5-polarization is expected to be very 
small, which makes the detection of 5-polarization quite difficult. In addition, the signal of RGWs 
in 5-polarization can be contaminated by the E-B mixing due to the partial sky coverage \\\ . 
beam asymmetry [2] and cosmic lensing effect \}\ . These all can degrade the detection ability of 
the "55" method. 



B. U CT case 



we have detailed discussed the 



Different from the "55" method, in the previous work \. 
" TE' method, detecting the signal of RGWs by the CMB TE power spectrum. In this method, 
the amplitude of Cj E is two order larger than Cf B . Another advantage of this method is that, the 
E-B mixing, which can occur for some reasons, nearly cannot influence TE power spectrum. So 
it cannot degrade of the detection ability of this method. However, in the previous work j^jj], we 
find that, the uncertainty of the TE estimator Dj E is very large, due to the cosmic uncertainty. 
So comparing with "55" method, "T5" method has not only the larger signal, but also the larger 
uncertainty. 

In this paper, we shall develop the "T5" method by combining the CMB TE and TT power 
spectra. In the real observations, the amplitude of Cj T is much larger than that of the other three 
power spectra. So combining the TE and TT power spectra are expected to be a more effective 
way to detect RGWs. We denote it as "CT" case. In this case, the likelihood function is 

Cc^X\f(Dj E ,DD. (18) 
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Using the pdf in (|13p . this likelihood function can be rewritten as 
- = + ^ - + >n(<l - ,?)(^""(^)} + C,( 19 ) 

C. "CTE" case 

In this case, in addition to the TE and TT power spectra, we shall include the i?-polarization 
power spectrum Cf E . By comparing with "CT" case, we can investigate the contribution of 
-E-polarization for the detection of RGWs. In this case, the likelihood function is 

C C TEOcl[f(Dj E ,Dj T ,Df E ). (20) 

i 

Using the pdf in (jlip . this likelihood can be written as 
" ' >■ An. = £ { ^ ( ^ + jfa ~ H) + -eff b (4(1 - Pf)(^-f) 2 ) } + <%W 

D. "CTEB" case 

This case will use all the CMB power spectra, CT E , Cf T , Cf E and Cf B , so it is a combination 
of "CTE" and "B". By investigating this case, we can determine the best constraint of RGWs by 
the CMB observation. In this case, the likelihood is 

£cteb oc J] f(Df E , Dj T , D EE , Df B ), (22) 

i 

which is the product of £cte and Cb, i-e. 

-£-CTEB = ^CTE-Cb- (23) 



IV. THE SIMULATION METHOD AND THE RESULTS 

As the previous work in this section, we shall use the maximum likelihood analysis, based 
on the simulated data, to discuss the sensitivities for the detection of RGWs in these four cases. 

Before proceeding on the simulation method, we shall firstly introduce the background cosmo- 
logical model. Throughout this paper, we shall adopt a set of typical cosmological parameters as 
follows 



h = 0.732, n b h 2 = 0.02229, Q m h 2 = 0.1277, Q k = 0, T reion = 0.089. (24) 
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Since in this paper, we focus on the detection abilities for the RGWs, in Sections IIVI and PVT we shall 
only consider the constraint on the parameter r. In Section |VH we shall extend to the constraint 
on the other three parameters nt, A s , n s , and discuss their influence on the constraint of r. The 
extent of the constraints on the cosmological parameters (h, f^, O m , Q^, r re j on ) remains an open 
question in this paper. Actually, by the forthcoming observation of Planck satellite, the constraints 
on these cosmological parameters are expected to be very tight. For example, the constraint on 
Q^h 2 would be Af^/i 2 = 0.00017, the constraint on r re j ori would be Ar re j on = 0.005 18|. which are 



expected to have little influence on the determination of RGWs. In all this paper, we take specific 
values for cosmological parameters as in (|24p and assume that they are perfectly known. 

A. The method 

In this section, we shall use the maximum likelihood analysis to investigate the constraint on 
the cosmological parameters. This method has been used in the previous work 29j] for the CMB 



analysis and in the work [40j] for dark energy analysis. If we consider the "B" case, the steps of the 
method can be listed as the follows (the similar steps can also be used in the "CT", "CT.E" and 
"CTEE 1 cases): 

Step 1 We build the pdf of the estimator: f(D BB ), which have been given in Eq. (|10p . 
Step 2 According to this pdf, we generate N sets of random samples {D BB \£ = 2,3, • • ■,£ ma , x } 
(we call each sample as a "realization"), where the input model has the parameters (r, fit, A s , n s ) 



46]. 



Step 3 We separate the these parameters into two sets: the first set includes the so-called 
unfixed parameters, and the second set includes the fixed parameters. In this section, we are only 
interested in the constraint on the amplitude of RGWs, so we consider the simplest case, where 
the only unfixed parameter is the tensor-to-scalar ratio r. The other three parameters, nt, A s and 
n s , are all the fixed parameters. In Section IVT1 we shall discuss the influence of other parameters 
on constraint of r, so we shall choose more than one parameters as the unfixed parameters. 

Step 4 We fix the fixed parameters as their input values and set the unfixed ones as the free 
parameters. Using Eq. (|17p. an automated search, which uses the numerical technique of simulated 



annealing 



411 ] . finds the maximum of likelihood C-q for each realization. 



Step 5 To measure the certainty with which the unfixed parameters can be determined, we 
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examine the distribution of the maxima from the simulations. 

Evaluation of the likelihood function on a fine grids of the unfixed parameters shows that the 
maximum found by the automated procedure differs negligibly from the true maximum. Performing 
N realization allows us to determine the standard deviations of the unfixed parameters with a 
fractional error of (2iV)- ' 5 . When N = 300, the fractional error is 4%, and when N = 1000, the 
fractional error is 2%. 



B. Results 

We apply the simulation method to the "B", "CT', "CTT" and "CTEH' cases. The values of 
l max , N and the input values of parameters are adopted as the follows: 

Imax = 100, N = 1000, r = 0.3, n t = 0.0, A s = 2.3 x 10~ 9 , n s = 1.0. (25) 

Since when considering the Planck instrumental noises, the contributions of RGWs in the CMB 
power spectra Cf x ' are important only at the large scale (£ < 100) 23], we have adopt £ ma , x = 100, 
i.e. only using the simulated data in the range £ < 100 in the likelihood analysis. For each case, 
we generate (N =)1000 realization, which makes that we can determine the standard deviation of 
the unfixed parameter with a fractional error of 2%. 

In this section, we only set r as the unfixed parameter, i.e. the spectral indices nt, n s and the 
amplitude A s will be fixed as their input values in the calculation. As mentioned above, for any 
two pivot wavenumbers ko and k%, the constraints on r(ko) and r{k\) are exactly same. So in the 
discussion in this section, we shall not differentiate r at the different pivot wavenumbers. 

FIG.l presents the distribution of the maxima r p in each case. We find in all of these realization, 
the values of r p are close to the input value f = 0.3. In "5" case, we find the average value of 
r p is rp~ = 0.305. In "CT' case, we have = 0.298. In "CTT" case, we have = 0.303, and 
in "CTEB" case, we have = 0.304. These four average values are all equal to the input value 
within 1.7% (< 2%) and hence there is no evidence for bias. 

In the different cases, the diffusion of r p is different. The standard deviation of r p in "5" 
case is Ar p = 0.067, so we can conclude that an experiment of this type can determine r with 
Ar = 0.067 ± 0.001. In "CT" case, the standard deviation is Ar p = 0.078, which is 16.4% larger 
than that in "5" case. So we conclude that "CT' is a little less sensitive than "5" for the constraint 
of r. In "CTE" case, the standard deviation is Ar p = 0.070, which is 10.3% smaller than that in 
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FIG. 1: The distribution of r p in the 1000 realization. The black dots (upper left panel) denote 
the result in "5" case, the red dots (upper right panel) are the result in "CT" case, the magenta 
dots (lower left panel) are the result in "CTTi" case, and the blue dots (lower right panel) denote 
the result in "CTEB 1 case. In all these panels, we have considered one free parameter r in the 
likelihood analysis. The input simulated data are up to £ max = 100, and the input value is r = 0.3. 

"CT" case. So including the E'-polarization data, the constraint on r becomes tighter. However, it 
is also 3.0% larger than the value of Ar p in "TT case. In "CTEW case, we have Ar p = 0.047. which 
is much smaller than the others. So we conclude that, taking into account all the simulated data, 
the constraint on the parameter r can be much improved. Comparing with "T?", the uncertainty 
of r is reduced by 29.8%, and comparing with "CT", the uncertainty is reduced by 39.7%. 

Similar to the discussion in our previous work 
for the RGWs, we define the signal-to-noise ratio 

S/N = r/Ar p , (26) 

where f is the input value of r. We can determine this quantity with different input r. For each 
input, we generate 1000 realization, and calculate the quantities r^, Ar p and S/N. The value of 
S/N as a function of f are shown in FIG. 2. From this figure, let us firstly investigate the detection 
abilities in the four cases. We find, in "5" case, the parameter r can be determined at 2a level 
when f > 0.06. In "CT" case, r can be determined at 2a level when f > 0.16. In "CT£" case, r 



23J, in order to describe the detection abilities 
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input tensor-to-scalar ratio: r 



1 



FIG. 2: The value of S/N depends on the input value r. The black, red, magenta and blue dots 
(lines) are the simulation (analytic approximation) results in "B" case, "CT" case, "CTT?' case 
and "CTEB 1 case, respectively. 

can be determined at 2a level when f > 0.13, and in" CTEB" case, r can be determined at 2a level 
when f > 0.05. 

From FIG. 2, We can also compare the sensitivities in the different cases by the values of S/N. 
Comparing "CT' and "5", we find the former one is more sensitive when f is large, and the latter 
one is more sensitive when r is small. "CT' is more sensitive than "5" when f > 0.5. "CTT7' is 
more sensitive than "B" , when f > 0.3. In order to investigate the contribution of T?-polarization 
on the detection of RGWs, we can compare the values of S/N in "CT' and "CTTi" cases. From 
the FIG. 2, we find the quantity S/N in "CTE' case is always 15% larger than that in "CT' case. 
So considering the TJ-polarization, the constraint on r can be improved for any f. From the FIG. 2, 
we also find that, as the combination of "CTE 1 and "B\ " CTEB" is more sensitive than the other 
three cases. When r is small, the sensitivity in "CTEB" case is close to that in "i?" case, since 
in this case, the sensitivity of "CTTi" is very weak. When f is large, the sensitivity in "CTEB" is 
close to that in "CTT?' case. 

We have also applied the simulation method to another condition: the input quantities are all 
exactly same with those in Eq. ([25]) . except the value of £ max - In this case we adopt £ max = 500, i.e. 
the simulated data Df x ' in the range t < 500 are used for the likelihood analysis. In Table [J we 
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TABLE I: The mean values and standard deviations of r p . In the likelihood analysis, we have 
considered one free parameter r. 



input £ max 


output parameter 


B 


CT 


CTE 


CTEB 


100 


rp~ ± Ar p 


0.305 ± 0.067 


0.298 ± 0.078 


0.303 ± 0.070 


0.304 ± 0.047 


500 


f^± Ar p 


0.302 ± 0.067 


0.301 ± 0.080 


0.302 ± 0.066 


0.302 ± 0.047 
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FIG. 3: The same graph as in FIGj2j the only difference is that, in this figure we have assumed 
the realistic noise power spectra N xx are 4 times larger than the Planck instrumental noises. 

summarize the output values Tp~±Ar p in "CT", "CT£", "CTEB 1 cases, where f = 0.3 is used. 
We find that, the results in this condition is V6ry close to those in the condition with l m ax — 

100. 

This result testifies that, when considering the Planck instrumental noises, the contribution of 
RGWs in the CMB power spectra C xx are important only at the large scale (£ < 100). 

To this point we have assumed that the noise power spectra only come from the Planck in- 
strumental noise. However, synchrotron and bremsstrahlung radiation, thermal emission from cold 
dust, and unsolved extragalactic sources also contribute to the anisotropy and polarization of radi 
ation. These contaminations can enlarge the effective "noises" of CMB power spectra 



13, 



29, 



42j. 



In order to estimate the effect of these contaminations on the constraint of r, in this paper, we 
only simply assume the foreground will degrade the noise a*L(w) by a factor 2. Therefore, we take 
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into account the effect of foreground contaminations by simply increasing Nf x to 4iV/ A . In this 
case, by the exactly same steps as the previous discussion, we recalculate the signal-to-noise ratio 
S/N by the simulation method, where different input values r are considered. The quantity S/N 
as a function of f in four cases are shown in FIG. 3. Let us firstly discuss the results in "5" case. 
Comparing with the results in FIG. 2, we find the detection ability in "5" case is much decreased. 
Only when f > 0.25, the signal of RGWs can be detected in 2a level. However, in "CT' and 
"CTT?' cases, the results are similar with those in the previous condition with only instrumental 
noises (FIG. 2). Comparing the sensitivities in "CT' and "CTTi" cases, we find the difference is 
very small, which suggests that the contribution of T?-polarization for the detection of RGWs is 
negligible. In "CTEB" cases, when r > 0.12, the signal of RGWs can be detected in 2a level. 

In the previous work [23j], we found that WMAP5 TE data induces the best-fit model with 
r ~ 0.48. From FIGs. [2]and[3j we find that RGWs with r = 0.48 will be presented nearly at 9a 
level in "CTEB" case, when the Planck instrumental noises are considered. If the assumed realistic 
noises are considered, it will be presented at 6a level. These are all much better than the results 
in "T£" and "SB" methods [23]. 

V. ANALYTIC APPROXIMATION OF S/N 

In Section IIVI using the signal-to-noise ratio S/N calculated by the simulated data, we have 
investigated the detection abilities for the RGWs in four cases ( "5" , "CT" , "CTT" , "CT5B" cases). 
In order to better understand this signal-to-noise ratio and get an intuitive feel for the results in 
Section HVl in this section we shall give a simple analytic approximation of the signal-to- noise ratio. 
Similar to the discussion in Section IIVI in the analytic approximation, we will also be interested 
in a single unfixed parameter, the tensor-to-scalar ratio r. Other parameters (nt, A s and n s ) and 
background parameters (h, Of,, • • •) are all assumed to be exactly known. 

A. Analytic approximation of the likelihood functions 

In order to present the analytic expression of the signal-to-noise ratio, we need to express the 
likelihoods as the simple functions of variable r. We notice that the exact pdfs in Eqs. (|10p . 
(fTTj) . (fT5|) . (fT4"j) , are all very close to the Gaussian function, especially when I S> 1 (see for instant 
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23j, |35[). Based on the Gaussian approximation of these pdfs, the likelihood functions in Eqs. (|16|) . 
(EOD, (HHD, ([22]) can be simplified as ([A3]), (TATOl) . (lATTl) . (IA12T) (see Appendix [£] for the details), 
which can be rewritten in a unified form as follows: 

(27) 



21n£(r) = ^^ 



AD XX ' (f 



1=2 XX' 

In "5" case, we have XX' = BB; in "CT' case, we have XX' = TE,TT; in "CT£" case, 
XX' = TE, TT, EE; and in "CTEB" case, we have XX' = TE, TT, EE, BB. In each case, the 
likelihood function L depends on the variable r only by the power spectra C xx ' . 

In general, ignoring the possible contribution from the (vector) rotational perturbations, the 
CMB power spectra C xx ' can be presented as a sum of two contributions: density perturbations 
and gravitational waves: 

C xx '(r) = C xx '(dp) + C xx '(gw), (28) 

where C xx (dp) and C xx (gw) are the contributions of density perturbations and gravitational 
waves, respectively. We should remember Cf B (dp) = 0. In the likelihood analysis, we have fixed 
the parameters rif, A s , n s as their input values, and only considered a single free parameter r. 
C¥ x ' (gw) depends on the variable r, which can be written as 

Cf X '(9w) = ^)c xx '(9w), (29) 

where C xx ' (gw) are the power spectra C xx ' (gw) at r = f. Inserting (|29p in (|28p . we get 

C xx '(r) = C xx '(d P ) + Q C xx '(gw). (30) 

Now, let us return to the likelihood function. Inserting (|30p into Eq. (|27p . we obtain that 

-2lnC(r) = J2J2»f X '^-^ X ') , (31) 

e=2 xx' 

where the quantities ol xx> and P xx ' are defined by 

n .xx>-( c? X '(9») Y oxx' _ Df x ' (f) - C xx ' (dp) 



After straight forward manipulations, the expression (|3ip can be rewritten as the following form 
- 2C(r 



l ">» / „ V^max V- n XX ' R XX ' \ 2 

(ES°f X, > [h T'~T a ** +C ' <33) 
1=2 xx' \ 2-^i=2 l^xx' a e J 

where the separate part C is independent of the variable r. In the following, based on this formula, 
we shall discuss the signal-to- noise ratio S/N . 
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B. Analytic approximation of S/N 



Now, let us investigate the likelihood in (|33p . First, we shall discuss the peak of the likelihood 
function. We notice that this likelihood is a "Gaussian form function" of the variable r. The 



maximum is at r p , which is 



XX' oXX' 



V'™r rv xx ft 



V'"- V ,,.YV " r * ^ 
2^,1=2 l^XX' a l 

From the definition of a xx ' and @/F x ' in Eq. (|32p . we find the value of r p not only depends on the 
values of Cf x '(f), the input model, and N xx , the noise power spectra, but also depends on the 
values of D¥ x , the simulated data. So for any two realization, even if generated by the exactly 
same input model and noises, they have the different values of r p . 

From the likelihood in Eq. (|33h . we can also obtain the spread of the likelihood Ar, which is 



Ar = — TV! ■ ( 35 ) 

XX' \ ' 



V^maxV^ XX ' 

2^=2 /—/XX' a £ J 

The spread of the likelihood only depends on the values of C? x (r), the input model, and N xx , 
the noise power spectra. So we get the conclusion, for any two realization, as long as they have 
the same input model and noises, they have the same Ar, the spread of the likelihood function. 
In the previous work j^, we have defined the signal-to-noise ratio S/N as 

S/N = f/Ar. (36) 

Note that, in order to distinguish from S/N defined in (|26p . we denote the signal-to-noise ratio 
in our previous work as S/N. In the following, we will find these two definitions have the exactly 
same values. Using the formula in ([35]) and the definition of af x ' in ([32]) . we obtain 



S/N 



This is the finial analytic result of the quantity S/N, which depends on the input power spectra 



431 ] . by the Fisher Matrix analysis, the 



Cf x (f) and the noise power spectra N xx ' . In the work 
authors have obtained a same result as Eq. (|37p in the case of XX' = BB (the result in "5" 
case). The results in Eqs. (|33]) , ([35]) and ([37]) describe the constraint on the tensor-to-scalar ratio 
r, based on one set of simulated data {D xx \l = 2, 3, • • •, 4nax}- 

However, in our discussion in Section IIV] we have considered another case. In the simulation 
method, based on a same input cosmological model and noises, we have randomly generated N 
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(N ^ 1) realization. For each realization, we can obtain a maximum of likelihood r p . From these 
r p , we have calculated the mean value and standard deviation of r p . From the simulation, we 
find the mean value Tp~ is close to the input value r, and the standard deviation Ar p stands for 
the uncertainty of the parameter r, in the likelihood analysis. Now, we shall prove these in the 
analytic approximation. 

In the analytic approximation, from Eq. (|34p we can also calculate the values of and Ar p . 
The mean value is 

r p ~ V ^ xx ' 1 j 

Considering the definition of (3f x in Eq. ()32p . the relation {D xx (f)) = C xx (f) and the Eq. 
(j30|) . we can obtain that {f3 xx ) = 1. Inserting it into Eq. ([38]) . we get a relation 

Tp = f. (39) 

This relation suggests that, the mean value of r p is equal to the input value r, which is consistent 
with the simulation result in Section IIV1 



We can also discuss the standard deviation of r p , which is defined by Ar p = y f p — Tp~ . Using 
Eq. (|34p and the relation (D xx ' (f)) = C xx> (f), after a straight forward manipulations, we obtain 
that 

Ar p = \T/2' ( 4 °) 

Comparing (|40p with (|35p . we find Ar p = Ar. By the formula (|40p . we can discuss the signal-to- 
noise ratio S/N, defined by Eq. (|26p . Taking into account the definition of af x in Eq. (I32p . we 
can write the signal-to-noise ratio as 



S/N = S/N 



2 



which depends on the input power spectra C xx (r) and the noise power spectra N xx . In Eq. 
(JUD, we should remember that, XX' = BB in "5" case, XX' = TE,TT in "CT' case, XX' = 
TE, TT, EE in "CT£" case and XX' = TE, TT, EE, BB in "CTEB' case. From the expression in 
(|4ip . we find that, the two definitions of signal-to-noise ratio, S/N and S/iV have the same values. 
They all stand for the detection abilities for the RGWs. 

Using Eq. (I4ip . and taking into account the corresponding noises power spectra N xx , in FIGJ2] 
and [3 we have plotted the quantities S/N as a function of f in "5", "CT', U CTE\ "CTEB 1 cases. 
We find they are all exactly consistent with the simulation results. 
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C. Understanding the analytic approximation of S/N 



Now, let us investigate the approximation formula (|4ip . which can be rewritten as 



(W . (W jgg(iPi) 2 (42) 

In this expression, C xx ' (gw) is the contribution of RGWs to the total power spectra, which 
determines the strength of the signal of RGWs. AD XX ' is the uncertainty of the estimator, which 
serves as the corresponding 'noises'. So the total signal-to-noise ratio S/N (or S/N) is determined 
by the sum of the ratios between RGWs signal and the corresponding 'noise' at every multipole 
and XX'. This result is consistent with that in method, which has been obtained in our 

previous work 23]. 



From Eq. (I42p . we can discuss the contribution of each I to the total signal-to-noise ratio. We 
define the signal-to-noise ratio at the individual multipole £, (S/N)i as below: 



XX' 



Thus the total signal-to-noise ratio can be written as the following sum: 

(S/N) 2 = (S/N) 2 = ^(S/N) 2 . (44) 



Let us discuss the quantity (S/N) 2 . Taking into account the noise power spectra N xx , and 
adopt the input value f = 0.3, in FIG. H] we plot the quantity (S/N) 2 , as a function of I. Let us 
firstly focus on the lines in "B" case (black lines). We find that (S/N) 2 is sharply peaked at I < 10. 



As mentioned in our previous work [23(], the main contribution in "B" case comes from the signal 
in the range £ < 10. Especially, by comparing the line in left panel (Planck instrumental noises are 
considered) with the right panel (noise power spectra are 4 times larger than Planck instrumental 
noises), we find when increasing the noise power spectra, the value of (S/N) 2 reduces by a factor 
2, at the scale i < 10. However, at the scale i > 10, the value of (S/N) 2 reduces by a factor 10. So 
we get the conclusion, when increasing the noise power spectra, in "5" case, the contribution from 
£ < 10 becomes more and more dominant. Since in the range t < 10, the BB power spectra Cf B 
are mainly generated by the cosmic reioniztion 8, j|, the sensitivity in "i?" case strongly depends 
on the cosmic reionization process. 

Let us turn to the the quantity (S/N)} in "CT' and "CT£" cases, which are plotted in red 
(dashed) and magenta (solid) lines in FIG HI We find that, in both panels of FIGHJ the quantities 
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Multipole: I Multipole: I 



FIG. 4: This figure shows the individual signal-to-noise ratio (S/N)f as a function of £. The 
presented in left panel is the result with the Planck instrumental noises, and the presented in right 
panel is the result in the case where we assume the realistic noise power spectra N xx are 4 times 
larger than the Planck instrumental noises. 

{S/N)} in "CT" and "CT£" cases, are all peaked at £ ~ 30. Among the total S/N, the main 
contribution comes from the intermedial range 10 < £ < 70. In the very large scale £ < 10 and the 
small scale £ > 70, the values of {S/N)} are all very small. Their contributions to the total S/N are 
negligible. Since {S/N)} is very small in the scale £ < 10, the influence of cosmic reionization on the 
detection abilities are also not obvious, which is same with the " TE' method, but different from the 
ll BE" method. By comparing the solid and dashed lines in the left panel with the corresponding 
lines in the right panel, we find that, increasing the noises, the values of {S/N)} in "CT' and 
li CTE' cases have no obvious change, which induces that the total S/N in "CT" and "CTi?' cases 
have no obvious change (see FIG [2] and [3j). 

From Eq. (@3J, we find, in "CTE 7 (or "CT") case, the quantity {S/N)} is a simple sum of 
the portions with XX' = TE, TT and EE. So we can also discuss their contributions to the 
quantity {S/N)} separately. In FIG03 we plot these three portions of {S/N)} in dash-dotted 
lines, dashed lines, and dotted lines. In this figure, the solid lines denote the sum of these three 
portions, which exactly corresponds to the magenta lines in FIGH1 We find, when considering the 
Planck instrumental noises, these three portions are close to each. In the range 10 < £ < 70, the 
largest contribution comes from the component XX' = TT, which is 2 or 3 times larger than the 
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Multipole: I 




FIG. 5: This figure shows three quantities of (C xx ' (gw)/ADf x ') 2 (XX' = TE, TT, EE) as a 
function of t. As the combination of these three components, the individual signal-to-noise ratio 
(S/N)} in "CTTi" case are plotted in solid lines. In the left panel, we have considered the Planck 
instrumental noises, and in the right panel, we have assumed the realistic noise power spectra are 
4 time larger than the Planck instrumental noises. 

components XX' = TE and XX' = EE. However, in the case with large noises (right panel in 
FIGJS]), the contributions of the component XX' = EE rapidly decreases, and becomes negligible 
among the quantity (S/N)f, which induces that the value of the quantity (S/N)f in "CT£" case 
are very close to that in "CT" case (see the right panel of FIG|4|). So the total S/N in "CTE" and 
"CT' cases are also very close to each other (see FIG|3]). 

We can also discuss the quantity (S/N)j in "CTEB 1 case, which are plotted in blue (dash- 
dotted) lines in FIGH From Eq. flSD, we find the quantity (S/N)j in "CTEB" case is a simple 
sum of (S/N)} in "CT£" case and in "5" case. In the range of I < 10, the quantity (S/N)} in 
"CTEB" case is close to that in ll B" case, and in the range of t > 10, it is close to that in ll CTE J 
case. 

VI. EFFECTS OF FREE PARAMETERS: n t AND A s , n s 

In the previous sections, by both the simulation and the analytic approximation, we have 
discussed the constraints of the tensor-to-scalar ratio r in "B", "CT", Ll CTE\ "CTEB" cases. 
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However, in the real detection, we always have to constrain all the cosmic parameters, including 
h the Hubble parameter, the baryon density, Q m the matter density, Qk the spatial curvature, 
T reion the reionization optical depth. The parameters also include the scalar spectrum parameters: 
A s the amplitude of scalar spectrum and n s the scalar spectral index, and the tensor spectrum 
parameters: r the tensor-to-scalar ratio and nt the tensor spectral index. 

As mentioned in Section IIVI in all this paper, we shall not consider the constraints of the 
background cosmological parameters, and assume they have been exactly determined. In the 
likelihood analysis in Sections IIVI and IVl we have considered the case with only one free parameter 
r. The other parameters nt, A s and n s are all fixed as their input values. Based on this assumption, 
we have discussed the constraint on r in "5", "CT", U CTE\ "CTEB" cases. Thus a question arises, 
if the parameters, nt, A s , n s are also set free in the likelihood analysis, whether they can influence 
the constraint on the parameter r. 

In this section, by the simulation method, we shall answer this question. In the Section IVI A| 
we shall discuss the constraints on the free parameters, r and nt, and investigate the effect of nt 
on the constraint of r. In Section TVI Bl we shall consider the free parameters r, nt, A s and n s , and 
investigate the effects of A s and n s on the constraint of r. 

A. Effect of the free parameter n t 

Here, we shall use the simulation method described in Section HV A[ We choose the parameters 
r and nt as the unfixed parameters, A s and n s as the fixed parameters. The values of £ ma . x , N and 
the input values of parameters are adopted as the follows: 

4iax = 100, N = 300, r = 0.3, h t = 0.0, A s = 2.3 x 10~ 9 , n s = 1.0. (45) 

The background cosmological parameters are adopt as in Eq. (|24p . We consider the Planck 
instrumental noises and Planck window function, which are given in Eqs. (|6ati6dp . N = 300 
suggests that our following simulation results (r^, Ar p , ni, Ant) have 4% statistical error. 

As mentioned above (n s = 1 and nt = 0), for any two different pivot wavenumbers ko and k\, 
the tensor-to-scalar ratio r{ko) and r{k\) have the different constraints, due to the free parameter 
nt- Although they have the same input values r = 0.3, due to the formula in ((3|) and fit — n s + 1 = 0. 
As the first step, in the likelihood analysis, we choose the pivot wavenumber 

k = 0.05Mpc _1 . (46) 
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FIG. 6: The locations of the maxima from 300 realization projected into nt — r p (upper panels), 
and nt — r* (lower panels) planes. Black (red, magenta, blue) dots denote the results in "5" ("CT", 
"CTE' , "CTEB") case. In all these graphs, we have considered two free parameters (r, i%t) in the 
likelihood analysis. The input simulated data are up to £ max = 100, and the sign " + " denotes the 
input values of the parameters. 

Presented in FIGJ6] (upper panels) shows the maxima projected into nt — r p plane from 300 real- 
ization. First, we discuss the "5" case. The result is Tp~ ± Ar p = 0.351 ± 0.235. The uncertainty 
of r becomes nearly four times larger than the previous one (the result in the case with fixed nt), 
due to the free tensor spectral index m- The constraint on n% is: nt ± Ant = —0.022 ± 0.240. 

From FIGEJ we also find the strong correlation between nt and r p , which can be easily un- 
derstood. It is due to we have chosen the pivot scale ko = 0.05Mpc _1 . However, the quantity r 
in this scale, r(ko), is not the one which is measured most precisely. We assume that there is a 
tensor-to-scalar ratio r*(k%) (the tensor-to-scalar ratio at the pivot wavenumber k%), which can be 
measured most precisely. We expect this quantity r*(&!) has no correlation with nt- In this paper, 
we call as the 'best pivot wavenumber'. Following Eq. ([3]), we can relate r*(kt) and r(ko) by 
the following formula 



Since in the calculation, we have adopted the input tensor spectral index fit = 0, we have r*{kl) = 




(47) 
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TABLE II: The best pivot wavenumber k+, the mean values and the standard deviations of (r, 
nt, r*). In the likelihood analysis, we have considered two free parameters (r, nt). 



input £ max 


output parameter 


B 


CT 


CTE 


CTEB 


100 


/^(Mpcr 1 ) 


1.26 x 10~ 3 


3.13 x 10" 3 


3.43 x 10~ 3 


2.25 x 10~ 3 


100 


r* ± Ar* 

1 p -J- L-i.1 p 


0.296 ± 0.072 


0.288 ± 0.081 


0.293 ± 0.073 


0.300 ± 0.049 


100 


W±An t 


-0.022 ± 0.240 


0.051 ± 0.475 


0.023 ± 0.369 


-0.003 ±0.166 


100 


r^±Ar p 


0.351 ±0.235 


0.775 ± 1.329 


0.502 ± 0.678 


0.335 ±0.167 


500 


^(Mpc- 1 ) 


1.52 x 10~ 3 


2.97 x 10~ 3 


3.75 x 10" 3 


2.37 x 10" 3 


500 


T;±Ar; 


0.291 ± 0.066 


0.291 ± 0.083 


0.293 ± 0.067 


0.298 ± 0.046 


500 


W t ±An t 


-0.016 ±0.226 


-0.024 ±0.433 


-0.019 ±0.327 


-0.009 ±0.159 


500 


Tp-±Ar p 


0.352 ± 0.249 


0.559 ± 0.950 


0.399 ± 0.414 


0.325 ±0.172 



f(fco) = 0-3. However, the uncertainties of these two quantities (r* and r) are expected to be 
different. 

We use the following steps to search for the best pivot wavenumber k% : 
Step 1 Randomly choose a pivot wavenumber k' , which is different from ko- 
Step 2 Calculate the value of r p (k') by the formula in Eq. ([3|). 

Step 3 Project the maxima of the likelihood functions for 300 realization into nt — r p (k') plane. 

Step 4 In nt — r p (k') plane, if r p (k') correlates with nt, we iterate the same steps from Step 
1. Otherwise, if r p (k') has the weakest correlation with n t , we get the result: k! = k*, and 
r p (k')=r* p (k* t ). 

By these four steps, we find, in "5" case, the best pivot wavenumber is k* = 1.26 x 10~ 3 Mpc _1 . 
In FIGH3 we plot the distribution of nt — r* (left lower panel). As expected, the correlation between 
nt and r* disappears. We also calculate the average value and standard deviation of r*, which is 
r* ± Ar* = 0.296 ± 0.072. The average value of r*, is equal to the input value r* = 0.3 within 
1% and hence there is no evidence for bias. The standard deviation of r* (Ar* = 0.072) is much 
smaller than that of r p (Ar p = 0.235), but close to the result Ar p = 0.067 gotten in Section ITVT 
where only free parameter r is considered. Hence we conclude that, if we choose the best pivot 
wavenumber, the free parameter nt cannot influence the constraint on the tensor-to-scalar ratio. 
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We can also consider the constraints on r and nt in the other cases. The distributions of r p 
and rit in the 300 realization are all plotted in FIG. 6 (upper panels). The strong correlations exist 
in all these panels. By the exactly same steps, we can find the best wavenumber k%, which are 
all listed in Table HI For example, in "CT' case fc t * = 3.13 x lO^Mpc" 1 and in U CTE' case 
k* t = 3.43 x H^Mpc -1 . In these two cases, the best pivot wavenumbers are close to each other, 
which are all much larger than that in "5" case. In "CTEB" case, the best pivot wavenumber is 
k* t = 2.25 x 10 -3 Mpc -1 , which is between that in "5" case and that in "CT£" case. In all these 
three cases, the values of Ar* are all very close to that of Ar p gotten in Section IIVI when only 
free parameter r is considered. Hence, we obtain the same conclusion, if we adopt the best pivot 
wavenumber, the free parameter m cannot expand the constraint on the tensor-to-scalar ratio (We 



should mention that, in the latter work |44|, we have completely proved this conclusion, and given 
the analytic formulae for the best pivot wavenumber k% and the uncertainty Ant). 

From Table [III we a ^ so nn d that, the uncertainty of rit is always very large. For example in 
li CTEE" case, the constraint is Ant = 0.166, which is fairly loose for the determination of the 
physical model of the early universe. 

We have also considered another condition, where a broader range (^ max = 500) simulated data 
Df x ' are used for the likelihood analysis. The results are all listed in Table HH As expected, we 
find that in this condition, the values of h% Ar*, Ant are all close to those in the previous condition, 
where only simulated data in large scale (i max = 100) are considered. 



B. Effect of the free parameters A s , n s 

In this subsection, we shall extend the discussion in Section IVI Al to the more general case, 
where we consider four free parameters: (r, rij, A s , n s ). By the simulated data, we can investigate 
the effects of free parameters A s and n s on the constraint of r. The steps are exactly same 
with that in the Section [VI Al For the simplification, we shall use the parameter A' s , defined by 
A' s = A s /2.3 x 10~ 9 , instead of A s . 

We notice that, the power spectra Cf B only depends on Pt(k), but not on P s (k). Since Pt(k) 
is determined by the parameters At = {rA s ) and nt, in "5" case we cannot constrain the separate 
parameters: r, nt, A' s and n s . So in this subsection, we shall not discuss the "S" case. 

We firstly consider the condition, where the values of £ max , N, and the input values of parameters 
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are adopted as in (|45p . We adopt the pivot wavenumber as in (|46p . The likelihood functions peak 
at (r p ,nt,A' s ,n s ). FIG J7] presents the maxima projected into n± — r p plane from 300 realization. 
We find the strong correlation between nt and r p exists, which is because we have used the pivot 
wavenumber k$ = O.OSMpc" 1 . The outputs Yp~ ± Ar p , ni ± Ant, ^ ± AA' S , raj ± An s in "CT", 
"CTE" , "CTEB 1 cases are all listed in Table Hill We find, due to the uncertainties of A' s and n s , 
the values of Ar p and Ant ar e all larger than the corresponding results in Table [TTJ 

Similar to Section IVI Al we can discuss r*(k%), the tensor-to-scalar ratio at the best pivot 
wavenumber fe£. Following Eq. ([3]), we can relate r*(fe^ ) and r(A;o) by the following formula 

r(h) = r*(fc?) [pj . (48) 

Since in the calculation, we have adopted the input spectral index fit = and h s = 1, we find 
f*=f = 0.3. However, the uncertainties of r* and r are expected to be different. 

We search for the best pivot wavenumber k\ by the exactly same steps, listed in Section fVI Ai In 
"CT" case, the best pivot wavenumber is k\ = 3.11 x 10~ 3 Mpc~ 1 . Based on this pivot wavenumber, 
we find Ar* = 0.127. Comparing with the result of Ar* = 0.081, where only two free parameters 
r and n t are considered, the value of Ar* increases by 60%, due to the free parameters A' s and n s . 
Since we have only used the simulated in the large scale I < 100 in the likelihood analysis, the 
uncertainties of A' s and n s are fairly large (see Table IIIip . This makes the value Ar * is obviously 
increased. 

We have also considered the condition, where £ max = 500 is adopted. We find the constraints 
on A' s and n s become much smaller: AA' S = 0.005 and An s = 0.008, and the constraint on r* 
becomes Ar* = 0.091, i.e. the influence of A' s and n s on the constraint of r* becomes much smaller 
(increasing the value of Ar* only by 10%). 

In "CTi?' and ll CTEB" cases, we have also investigated the effects of free parameters A' s and n s 
on the constraint of r*. The results are all similar with those in "CT" case. The nt — r p and nt — r* 
planes are all plotted in FIGd We find in both cases, r p strongly correlates with m- However, 
as expected, r* does not correlate with nt- The best pivot wavenumber and the constraints 
of the parameters are all listed in Table IIIIl Based on these, we conclude that: In the likelihood 
analysis, if we only consider the simulated data in the large scale (I < 100), the constraints of 
r* and nt become much looser, due to the uncertainty of n s and A' s . However, if we considered 
the simulated data in the larger range {t < 500), the constraints on r* and n% only increase by 
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FIG. 7: The locations of the maxima from 300 realization projected into nt — r p (upper panels), 
and nt — r* (lower panels) planes. Red (magenta, blue) dots denote the results in "CT' { LL CTE\ 
"CTEtF) case. In all these graphs, we have considered four free parameters (r, nt, A' s , n s ) in the 
likelihood analysis. The input simulated data are up to £ max = 500, and the sign " + " denotes the 
input values of the parameters. 

~ 10%. Expectable, in the likelihood analysis, if the simulated data in the range I < 2000 (the 
real TT, TE and EE data, especially the TT data, in this range are expected to be well observed 
by Planck satellite 



181] ) are used, the influence of A' a and n s on the values of Ar* and Ant will 



become negligible. 



VII. CONCLUSION 

The upcoming observations of Planck satellite provide a very possible opportunity to detect 
RGWs in the CMB power spectra. In this paper, by both the simulation and the analytic approx- 
imation methods, we have discussed the detection abilities for RGWs in four ("5", "CT", "CTT", 
"CT£5") cases. The main conclusion can be summarized as: 1) In "5" ("CT", "CTT", "CTEB") 
case, the Planck satellite can detect the signal of RGWs at 2a level when r > 0.06 (r > 0.16, 
r > 0.13, r > 0.05). 2) Comparing "CTT" with "5", we find that, when r > 0.3, the value of the 
signal-to-noise ratio S/N is larger in "CTTi" case, and when r < 0.3, the value of S/N is larger 
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TABLE III: The best pivot wavenumber k%, the mean values and the standard deviations of (r p , 
n t , r* A' s , n s ). In the likelihood analysis, we have considered four free parameters (r, m, A' s , n s ). 



input £ max 


output parameter 


CT 


CTE 


CTEB 


100 


/^(Mpcr 1 ) 


3.11 x 10~ 3 


3.66 x 10- 3 


2.33 x 10~ 3 


100 


r* ± Ar* 


0.291 ±0.127 


0.295 ±0.115 


0.300 ± 0.069 


100 


ru±An t 


-0.003 ± 0.509 


0.009 ± 0.425 


-0.010 ±0.176 


100 


J>±Ar p 


0.783 ± 1.693 


0.546 ± 0.838 


0.331 ±0.178 


100 


nj ± An s 


1.006 ±0.052 


1.004 ± 0.050 


1.001 ±0.043 


100 


W S ±AA' S 


1.016 ±0.092 


1.012 ±0.090 


1.006 ±0.078 


500 


^(Mpc- 1 ) 


3.02 x 10" 3 


3.35 x 10~ 3 


2.22 x 10~ 3 


500 


r* ± Ar* 


0.290 ± 0.091 


0.294 ± 0.078 


0.297 ±0.051 


500 


nt± An t 


-0.024 ± 0.442 


-0.053 ±0.291 


-0.016 ±0.155 


500 


Tp-±Ar p 


0.570 ± 0.990 


0.341 ±0.341 


0.316 ±0.166 


500 


nj ± An s 


1.000 ± 0.008 


0.999 ± 0.008 


0.999 ± 0.007 


500 


aZ±aa' s 


0.999 ± 0.005 


0.999 ± 0.005 


0.999 ± 0.005 



in "5" case. If the realistic noise power spectra of Planck satellite is enlarged for some reasons, 
the value of S/N in "i?" case will be much reduced. However, in "CTi?" case, the value of S/N 
is little influenced. 3) The value of S/N is much larger in "CTEB J case than that in "5" case, 
especially when r > 0.1. 4) The free parameters nt, n s and A s , cannot reduce the value of S/N, if 
we consider the data in a large range and adopt the best pivot scale. 
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APPENDIX A: GAUSSIAN APPROXIMATIONS OF THE LIKELIHOOD 
FUNCTIONS 



In this appendix, by using the Gaussian approximation of the pdfs for the estimators D¥ x , we 
shall simplify the exact likelihood functions, given in Section IIIII 

1. Approximation of Cb 



First, let us focus on the analytic approximation of Cb- We use the following Gaussian function 



fc{D BB ) to approximate the exact pdf f(D BB ) 



fG{Df B {r)) 



1 



— exp 



^2ttAD bb (r j 

Inserting this formula into Eq. (|16p . we obtain that 

i r 

1 



{Df B {f)-C BB ? 
2{AD BB (r)) 2 



35] 



c B {r) = c n 



2irAD 



(D BB (f)-C BB ) 2 
2{AD BB {r)f 



(Al) 



(A2) 



where C is the constant for the normalization of the likelihood function, Df B (f) is the data, based 



on the input tensor-to-scalar ratio r. AD B (f) is the standard deviation of D 



*>BB, 



r). 



We should 



mention that, as a kind of approximation, Eq. (|A2|) can give results consistent with the exact 



likelihood function (the detailed discussion can be found in 



Up to a constant, we can rewritten the likelihood (|A2j) as follows: 



21n£ B (r) = £ 



fC BB 



1=2 



D BB (f) 



AD 



BB I 



(A3) 



which includes the variable r only by the power spectrum C 



BB 



2. Approximation of £cte 



Before proceeding on the "CT" case, let us firstly focus on the analytic approximation in "CTi?" 
case. The likelihood function Cqte depends on the pdf f(Dj E ,fJ T ,Cf E ), which is the Wishart 
function in Eq. (jlip . Similar to the approximation of f{Df B ), here we shall use the following 
multivariate normal function to approximate the exact Wishart distribution function: 



f G (D e ) 



(2^)3/ 2 |S| 1 A 



■ exp 



-hDe-difn-^Dt-dt) 



(A4) 
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where the vectors Dg and Ci are defined as Dg = [Dj E (f), Dj T '(f), Df E (f)] T , Q = 
[Cj E , Cj T , Cf E ] T . £ is the covariance matrix of the variable D^. Based on the Gaussian as- 



sumption of the CMB field, the estimators D xx ' have covariance as below (see for instant 
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covf^ D™) ~ 2{CJT + NJTW " ? (A5a) 

cm(D t ,D e ) - {2 £+l)f Bky ' {A5a) 

co V (D EE D EE ) - ^ EE + N EE W^f 

cov(D £ ,D e ) - {2i + 1)fsky > ( A5b ) 

covpr, - ( ^ )2 + (Cr + g^^ (A5c) 

-(^n ~ (ah) 

OT(Dp> D?T) = 2C, (|_±^ > , (A5e) 

C0V(D , £ Df£) _ zc, (|_^_ s ) (A5f) 

In order to investigate the cross relation between the estimators, we define the cross-correlation 
coefficient as 

PXX'YY' = j (AbJ 



/ cov( J Df x ', D xx ')coy(Dj DJY') 
From the relations in Eqs. (| A5all A5fj) . we can obtain that 



Pttee = Pi j Ptett = Pteee = Pe\ — — (A7) 

V 1 + Pi 



where pi is expressed in (|12p . which have been detailed discussed in our previous paper [23j. Taking 



23(]. 



into account the Planck instrumental noises, in the large scale (t < 100), we have pi < 0.45 
This makes that the correlation coefficients pttee, Ptett and pteee are all much smaller than 
1. So in the analytic approximation, we ignore the correlation between different estimators. Based 
on this approximation, we can simplify the multivariate normal function /g(-D) in (|A4j) as the 
following form: 

f G (D e )=Hf G (D? x '(f)), (A8) 

XX' 

where XX' = TE,TT, EE. The function fQ(D xx> (f)) is the following Gaussian function 

[Df x \r)-Cf x 'f 



2(l\D xx '(f)f 



(A9) 
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Inserting the approximation pdf (|A8P into Eq. (|20p and ignoring the independent constant, we 
get the approximation likelihood function, 



21n£ C TE(r) = £^ 



=2 XX' 



a 



XX' r)XX' 



Df x '(f) 



AD xx '(r) 



(A10) 



where XX' = TE, TT, EE. 



3. Approximation of £ct 



Let us turn our attention to the analytic approximation of Cqt- Similar to the discussion of 
£cte, we can get the approximation form of Cqt- Up to a constant, the likelihood is written as 



21n£ CT (r) = £]T 



=2 XX' 



c xx'_ D xx' { 
AD xx '{f) 



(All) 



where XX' = TE, TT. 



4. Approximation of £cteb 



We can also discuss the approximation form of likelihood £cteb- Since £cteb 
using Eqs. (IA3h and (IA10p . we obtain that 



- 21n£cTEB(^) 
where XX' = TE, TT, EE, BB. 



EE 

1=2 XX' 



c xx'_ D xx' 



if) 



AD XX '( 



^CTE-Cb, 



(A12) 
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